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Abstract. This is a survey-like paper presenting some of the recent combinatorial con- 
siderations on r-Dyck paths, r-Parking functions, and especially the r-Tamari lattices. 
Giving a better understanding of the combinatorial interplay between these objects has 
become important in view of their (conjectural) role in the description of the graded 
character of the S„-modules of bivariate and trivariate diagonal coinvariant spaces for 
the symmetric group. 



In this paper, our intent is to describe the interplay between various combinatorial con- 
structions that have recently been considered in relation with the study of "diagonal 
coinvariant spaces" for the symmetric group. This story started back in the early 1990's 
with the seminal work of Garsia and Haiman on the bivariate diagonal coinvariant space, 
and recently unfolded in a new direction with the work of the author on the corresponding 
trivariate case. Right at the beginning, 20 years ago, both r-Dyck paths and r-parking 
functions where involved in interesting ways in setting up formulas for the character of 
the relevant spaces, in the bivariate situation. Since then, many new developments have 
eloquently confirmed the role of these combinatorial object in this context. Still the saga 
of the study of these spaces is far from over, and one the main remaining open prob- 
lem is to prove the the "Shuffle Conjecture" of [11]. This conjecture gives an entirely 
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1. Introduction 
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combinatorial description of the graded character for the bivariate case, all in terms of 
r-parking functions. 

Very recently, it has become apparent (see [5]) that the r-Tamari order is a main player 
in the extension of this saga to the trivariate case. In particular, it allows for an exten- 
sion to the trivariate case of the Shuffle Conjecture. It is also worth mentioning that 
setting up this new conjecture has lead to interesting, and hard to prove, new combi- 
natorial identities (see [7, 9]). It seems interesting to present, in one unified text, the 
various combinatorial developments that play a role in this story, as well as the necessary 
background so the story be mostly self-contained. Several new research directions are 
discussed. 



2. Dyck paths and Tamari order 



Catalan numbers 



1 



2 72 



n+l\n 

1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 



(1) 



have a long an interesting history (see [19]) going back to Ming and Euler, in the mid 
18th century (a long time before Catalan). For Euler they appeared as the number of 





5 4 5 4 5 4 5 4 

Figure 1. Triangulations of the pentagon. 

triangulations of a polygon (see Fig. 1), and he gave both a formula, equivalent to (1), 
and the generating function 



1 - 2x- V / T^4 



x 



2x 2 



1 + 2x + 5x 2 + Ux 3 + 42x 4 + 132x 5 + 



(2) 



Still before Catalan, these numbers where generalized by Fuss ([12]) by considering other 
polygon decompositions. Following a growing tradition, we say that the resulting num- 
bers 

(ir + 1) 71 s 
n 



(j(r) ._ 



(3) 



r n + 1 

are the Fuss-Catalan numbers. There is a very extensive classical literature on a 
large number of combinatorial interpretations of these integers (see [19]); and many new 
interesting aspects are still being discovered. One especially interesting direction is a 
connection with the study of Shi arrangements, and their analogs for other Coxeter 
groups (see [1, 2]). 
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Observe that Cn^ may be considered as a degree n — 1 polynomial in the "parameter" 
r. Indeed, directly from Formula (3),we see that 

-(r n + 2)(rn + 3) • • • (r n + n). (4) 



so that 



nl 



a 
a 
c 
a 



(r 



(2r + 2)/2, 

(3r + 2)(3r + 3)/6, 

(4r + 2)(4r + 3)(4r + 4)/24, 

(5r + 2)(5r + 3)(5r + 4)(5r + 5)/120, 



At any r e N, these polynomial evaluate to positive integers. Thus, with r = 2, we get 
the sequence of integers 



1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715, 



The generating series 



x) 



±—( 

n=0 V 



(r + 1) n 



n 



x 



for Fuss-Catalan numbers, satisfies the algebraic equation 

Y(x) = 1 + xY(x) r+1 . 
Equivalently, the Fuss-Catalan numbers satisfy the recurrence 

E 

n=ni+n2+...+n r +i 
»(r) 



/-»(r) 



n{r)(~i{r) _ _ _ ^~<(r) 



(5) 

(6) 
(7) 



with initial condition C, 



(<■) 



1. 



We here consider Fuss-Catalan numbers as enumerators of r-Dyck paths. Recall that 
these are the paths in N x N, i.e.: sequences of "south" and "east" steps 1 , that 

• start at (0, n) and end at (rn, 0), and 

• stay below the line y = (—l/r)x + n, of slope — 1/r. 

For instance, we have the 2-Dyck path of Figure 2. More precisely, a r-Dyck path a is 
a sequence 

a = (po,Pi, • • • ,Pjv) 
of iV = rn + 1 points = (xi, yi), in N x N, such that 

• Po — (0, n) and p N = (r n, 0), 

• Pi + i = Pi + (—1,0) (a south step), or p i+1 = Pi + (0, 1) (an east step), and 
1 We chose this representation to simplify our later discussion. 
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(0,5) 




(10,0) 



Figure 2. The 2-Dyck path encodes as 00367. 



Xi < r (n — i/i), for all 1 < i < n. 



We say that n is the height of the path a, and we denote by 2?1 the set of r-Dyck 
paths of height n. It is well known that the number of r-Dyck paths of height n is given 
by the Fuss-Catalan number. 

For us, it will be practical to also describe Dyck paths in terms of the number of "cells" 
that lie to the left of south steps. More precisely, we bijectively encode height n Dyck 
paths in terms of sequences a^a^ ■ ■ ■ a n , with a, the ^-coordinate of the leftmost point on 
the path, among those that lie at height %. In this encoding, a = a\a^ ■ ■ ■ a n corresponds 
to a r-Dyck path if and only if 

(1) < cii < 0,2 < . . . < a n , and 

(2) for each i, we have < r(i — 1). 

For example, the 2-Dyck paths of height 3 are now encoded as 

000, 001, 002, 003, 004, 011, 012, 013, 014, 022, 023, 024. 

In such a description, each integer k, between and n (r — 1), occurs with some multi- 
plicity 2 nfc. Removing the 0-multiplicities, we obtain the (multiplicity) composition 

7(a) of the sequence a. Hence, we have 

7 (000) = 3, 7(001) = 21, 7(011) = 12, 7(002) = 21 
7(012) = 111, 7(022) = 12, 7(003) = 21, 7(013) = 111 
7(023) = 111, 7(004) = 21, 7(014) = 111, 7(024) = 111. 

By extension, we also say that 7(a) is the composition of the corresponding Dyck path. 
Clearly 7(a) is a composition of n (it sums up to n). 

The area, defined for ol in Dn ^ as 



area(a) := 



n 



i=i 



(8) 



2 possibly equal to 0. 
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is one of the well known interesting statistics on Dyck paths. Summing over the whole 

(r) 

set T>n , we obtain the g-enumerator polynomial 

C&to-E^' (9) 

a 

which satisfies the recurrence 

c2i(?) = E g r -^v r - 1)n2 ^ • • • cw, (io) 

n=ni +712 + . . . +n r+ 1 

with initial condition Cq = = 1. Equivalently, the generating series 

oo 

CW^^E^l^)^ (11) 
n=0 

for these g-Fuss Catalan numbers, satisfies the g-algebraic equation 

Y(x) = 1 + xY(q r x)Y(q r - 1 x) ■ ■ - Y(qx)Y(x). (12) 

For r = 1, we have 

C x (g) = 1 

C 2 (q) = 1 + q 

C 3 {q) = l + 2g + g 2 + g 3 

C 4 (q) = l + 3g + 3g 2 + 3g 3 + 2g 4 + g 5 + g 6 

C 5 (g) = l + 4g + 6g 2 + 7g 3 + 7g 4 + 5g 5 + 5g 6 + 3g 7 + 2g 8 + g 9 + g 10 




Figure 3. Dyck paths weighted by q' 



Tamari Order. The "r-Tamari poset" is obtained by considering the following order 
on r-Dyck paths. We say that ajflj+i . . . , Ofc is a primitive subsequence, of a r-Dyck 
path a\ . . . a n , if 

(1) a,j — a,i < r(j — i) for each % < j < k, and 

(2) either k = n, or a fc+1 — a, > r(/c + 1 — i). 



G 
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For each i, there is a unique such primitive subsequence. It corresponds to the portion 
of the r-Dyck path that starts at (cii,i — 1) and ends at the "first return" of the path 
to the line of slope 1/r passing through the point (a.j,z — 1). Whenever % is such that 
< eij, we set a < (3, where 

(3 :=ai aj_i (aj - 1) ... (a k - 1) a k+1 . . . a n , 

with aj . . . a k equal to the primitive subsequence starting at a,. Recall from [5] that the 
r-Tamari order 3 is the reflexive and transitive closure of this covering relation a < 
Its largest element (denoted, in the usual manner, by 1) is the Dyck path encodes as 
00 • • ■ 0, and its smallest (denoted by 0) is the Dyck path for which Oj = r (i — 1). For 
n = 3, and r = 4, we get the poset of Figure 4; and, for n = 4 and r = 1, the poset of 



008 




044 



048 



Figure 4. The 4-Tamari poset, for n = 3. 

Figure 5. 

It has recently been shown (see [22]) that the r-Tamari poset is "EL-shellable" . This has 
a number of interesting topological consequences. We also have the geometric realization 

(2) 

of Figure 6 for the poset T>\ . This figures suggests that there is a natural polytopal 

(r) 

decomposition of T> n , with each polytope a Minkowski sum of polytopal realizations of 
1-Tamari posets (associahedron). 

3. Intervals 

As usual, in Poset Theory, an interval [a, /3] is the set of 7 € T>n^ such that a < 7 < (5. 
Thus, for [a, 0\ to be non-empty, we must have a. < (3. Recalling the classical recursive 



3 As far as we know this poset had not yet been considered, before [5]. 
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Figure 5. The 1-Tamari poset (associahedron), for n — 4. 




Figure 6. Looking up the 2-Tamari poset, for n — 4. 

definition of Mobius function on intervals, we have 

{1 if a = p, 

- V(a,x) otherwise. ( 13 ) 
a<x</3 

In particular, fi(a, 0) = if a ^ /3, since the sum is empty in that case. The EL- 
shellability of the r-Tamari poset implies that its Mobius function takes values in the 
set { — 1,0,1}. In [22], this Mobius function has been calculated for intervals having 
homotopy type of spheres, and these are explicitly characterized. 

Our aim now is to recall enumeration results for (non-empty) intervals in Vn\ For 
example, for n = 3 and r — 1, there are 13 intervals in the corresponding Tamari poset 
(which appears as the front face of Figure 5, with labels obtained by dropping the final 
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3): 

[012,012], [012,011], [012,002], [012,001], [012,000], 

[002,002], [002,001], [002,000], 

[011,011], [011,000], 

[001,001], [001,000], 

[000,000]. 

In [10], it is shown that, for r = 1, the number of intervals is 

2 An + T 



n (n + 1) \ n — 1 

Computer experiments led the author to conjecture more generally (see [5]) that X n (r), 

(r) 

the number of intervals in T> n , is given by the formula 

z»M = 4^U (r + 1)2 ? + r ), <«> 

n [rn + 1) \ n — 1 / 

with x{~ ) denoting the function that takes value 1 or 0, depending on whether or not its 
argument is true or false. Observe that X n {r) is a polynomial in the parameter r, indeed 

Mr) = 1, 

Mr) = (r + l)(r + 2)/2, 

Mr) = (r + l)(r + 2)(3r 2 + 7r + 3)/6, 

Z 4 (r) = (r + l)(r + 2)(4r 2 + 9r + 3)(4r 2 + 9r + 4)/24, 

Z 5 ( r ) = (r + l)(r + 2)(5r 2 + llr + 3)(5r 2 + llr + 4)(5r 2 + llr + 5)/120. 

Formula (14) has been shown to hold by Bousquet-Melou et al. in [7]. Their approach 
exploits the fact that the generating function for the sequence (14), may be described in 
terms of the series z{t) that satisfies the functional equation 

z{t)=t{l-z{t))-^ +2 \ (15) 

which expands as 

; W = Vl( r ( r + 2 >" + "-V (W) 
' n\ n — 1 / 

n>l v 7 

It seems that it would be natural to refine this enumeration of intervals to take account 
of their homotopy type. Once again the results of [22] strongly suggest that would get 
nice enumeration formulas in this manner. 



Length counting. For a < (5 in r-Tamari order, let us denote by d(a, ft) the length of 
the longest chain going from a to /3, and we say that d(a, ft) is the length of the interval 
[a,/3]. Observe that 

Lemma 1. For all a < ft, we have 

d(a,ft) < area(/3) — area(a). (17) 



Proof. We need only observe that this is trivially so in the case of the covering relation. 

□ 
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We may refine our interval counting, to take into account this notion of interval length, 
together with the area. Namely, we consider the polynomial 

a<(3 

Thus, we get 

V?{q,t) = 1, 
V$\q,t) = (1 + *)+?, 

V { z\q,t) = (l + t + 2t 2 + t 3 ) + (2 + t + 2t 2 )q+(l + t)q 2 + q 3 , 
V { ^(q,t) = (l + t + 2t 2 + 3t 3 + 3t 4 + 3t 5 + t 6 ) 

+ (3 + 2t + 4t 2 + 5t 3 + 4t 4 + 3t 5 )g 

+ (3 + 2t + 4t 2 + 3t 3 + 3t 4 )g 2 

+ (3 + 2t + 2t 2 + 3t 3 )g 3 

+ (2 + t + 2t 2 )q 4 + (l + t) q 5 + q 6 . 

Observe that this polynomial is not symmetric in q and t, but we clearly have 

DM( 9 ,0)=DM(0, g ). 

Another interesting experimental observation, about the distance statistic, is that we 
seem to have the identity 

aeV ( n r) agX>£ r) 



4. Parking functions 



For a sequence of positive integers p = a± a 2 ■ ■ ■ a n , let 

i{<p) := 6i & 2 ■■■b n (18) 

be the increasing rearrangement of ip. This is to say that hi < for all 1 < % < n. 
For example, for <p = 402040101, we have i(<p) = 000011244. 

One says that <p is a r-parking function if b\ ■ • ■ is such that bk < r(k — 1), for 
all k. We denote by Vn the set of r-parking functions. For example, we have the 
cardinality 49 set V% containing: 



000 


001 


002 


003 


004 


010 


011 


012 


013 


014 


020 


021 


022 


023 


024 


030 


031 


032 


040 


041 


042 


100 


101 


102 


103 


104 


110 


120 


130 


140 


200 


201 


202 


203 


204 


210 


220 


230 


240 


300 


301 


302 


310 


320 


400 


401 


402 


410 


420. 
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Perforce, i(<p) is a r-Dyck path. It is said to be the shape of ip. We denote by V a the 
set of all parking functions of shape a. 

In view of the above definition, we clearly have a S n -action (on positions) on the set Vn^ 
of r-parking functions. The sets V a decompose this action into disjoint orbits, and S n 
acts transitively on each V a . If 7(a) = C\02 ■ • ■ is the multiplicity composition of a 
r-dyck path a, then the Young subgroup 

is the fixator subgroup of a-shaped parking functions. It follows that the number of 
parking functions of shape a is given by the multinomial coefficient 

#v a = (") = — v (19) 

\j{oc)J c 1 \c 2 i---c k \ 

Parking functions may also be considered as coset representatives of the subgroup H := 
aZ, with u = (1, 1, . . . , 1), of the abelian group Z£ n+1 . More precisely, one shows that 

Lemma 2. Each coset of H, in Z™ n+1; contains a unique r-parking function. 

It follows that the set of r-parking functions of length n, denoted by Vn \ has cardinality 
{rn + l) n_1 . Thus, we get the identity 

(™+ir'= E t)) < 20 > 



More enumeration formulas. Motivated by an algebraic context discussed in sec- 
tion 6, and following computer experiments, the author was led to conjecture (see [5]) 
that the following elegant formula holds 

E x(«</3)( 7( ^ } ) =(r + ir(^ + ir 2 , (21) 

a,/3eX> ? ( l r) 

In other words, this is the cardinality of the set 

Q« := {(«, ip) I <p G P« and a < i(<p)}. (22) 

This was subsequently shown to be the case in [9] . There approach exploits the fact that 
the exponential generating function, for the sequence (21), may be described as 

J2(r + l) n {rn + l) n ~ 2 - = (1 - r z(t)) e (r+1)z{t \ (23) 

n>0 

where z(t) is the series 

v 7 n>0 

which satisfies the functional equation 

z (t) =te <r+l)z{t)_ (25) 
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Partitions and Young diagram. Recall that (integer) partitions A = A1A2 • ■ ■ 
of d, are sequences of integers Aj G N such that 

(1) Ai > A 2 > ■ ■ ■ > A fc > 0, and 

(2) Ai + A 2 + . . . + X k = d. 

We write A h d, when A is such a partition of d. Partitions of integers are often presented 
in the form of a Young diagram (see 1 Figure 7), for which we typically use the same 
notation. Recall that this is the following set of cells of N x N 

A = {(z,j) I 0<i<A i+1 }. 



1 1 1 1 1 1 1 

n 

FFhrn- 


1 




□□[□ 
□n B 

Eb @ 

il 


FFFEn 

Fj 

1 1 1 1 

z 


P 

~l 1 

1 1 1 1 

=FFh [ 
: ^ f 




lllll 

fln 

1 1 1 
ffl 


rfn 

m" u 1 

\zn 3 !□ 


to 




1 !:□ =3 :□ 






1 1 


: =1 ^ En 


□ EB Eb E 



Figure 7. Young diagrams of partitions of n < 8. 

It follows that we have an inclusion order on partitions, considering them as subsets 
of N x N. For two partitions A and fi, the skew diagram (or skew partition) A//i is 
simply the set difference 

X/fi := I E A and £ /i}. 

Clearly, we must have fi C A, for A//i to be non-empty. 

We may consider r-Dyck paths as the outer boundary of partitions A that are contained 
in the r-staircase shape 

#> = (r(n-l),r(n-2), ... ,r,0), 

together with the necessary segments along the axes so that we link (0, n) to (to, 0). 
Such a path is represented in Figure 8, with r = 1, n = 8, and A = 66531100. Here 
we add the necessary number of 0-parts to turn A into a n-part partition (contained in 
5^). This allows us to keep track of n, and makes easy the correspondence between such 

4 This figure is from Wikipedia's page on integer partitions 
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partitions and our previous encoding of r-Dyck paths. Indeed, this is but the increasing 
reordering of the parts of A. 
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Figure 8. The standard tableaux encoding of 35661010. 

Now, we denote by A + l n the partition that is obtained by adding one cell in each of 
the n rows (possibly empty) of A. Hence, the skew shape (A + l n )/A is a vertical strip, 
i.e.: with exactly one cell in each row. In formula, the corresponding set of cells is (see 
Figure 8) 

(A + l n )/A = {(z,A m ) | 0<i<n-l}. 

Exploiting this point of view on r-Dyck paths, shape A parking functions may simply be 
considered as standard tableaux 

r:(A + l")/A^{l,2,...,n}. 

Recall that this means that, for i < j, we must have 

r(i, A i+ i) < r(j, A i+i ), if A m = X j+1 . 

To get back our previous parking function description from this tableau description, we 
simply associate to each entry (of the tableau), the number of cells of A that lie (to the 
left of the entry) on the row in which this entry appears. This process may best be 
understood by considering an example such as that of Figure 8. 

Standardization. The classical process of standardization associates a permutation 
st(w) to each length n word w = v\ ■ ■ ■ v n , on an ordered alphabet. This permutation is 
obtained by successively replacing the letters in w by the numbers from 1 to n in the 
following manner. We first replace all instances of the smallest letter, say there are k of 
them, by the numbers from 1 to k going from left to right; then we proceed likewise with 
the next largest letter, replacing its instances by the numbers k + 1, k + 2, etc.; and so 
on until all letters have been replaced. For example, on the alphabet of the integers, 

st(97750032) = 86751243. 

Applying this standardization process to parking functions, we get a surjection 

st : § n , 

since permutations appear as special cases of parking functions (up to a shift). The 
partition of the set Vn \ whose blocks are the fibers of this surjection, has interesting 
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properties that seem to be tied to the Shi arrangement. Let us denote by Tr a ' the block 
associated to a given a in § n . To see why this notation convention may be considered as 
"natural", observe that the block of the identity corresponds to the set of r-Dyck paths. 

(2) 

2, the resulting partition of V% consists in the 6 blocks 

{000, 001, 002, 003, 004, 011, 012, 013, 014, 022, 023, 024}, 
{010, 020, 021, 030, 031, 032, 040, 041, 042} , 
{100, 200, 201, 300, 301, 302, 400, 401, 402} , 
{101, 102, 103, 104, 202, 203, 204} , 
{110, 120, 130, 140, 220, 230, 240} , 
{210,310,320,410,420}. 

Back to the r-Tamari order. In terms of this tableau description, we can extend the 
r-Tamari order to r-parking functions. Going up the order consists in moving to the left, 
and possibly merging, "correctly chosen" consecutive columns of r. 

The first rule characterizes the set of consecutive columns that may be moved in r of 
shape (A + 1™)/A. We first (freely) select a non-empty column of r, and draw the slope 
— 1/r line that stars a the top left corner of its lowest cell. Going right, this lines must 
necessarily hit back the r-Dyck path A (for the first time). We move one step to the left 
all columns between (and including) the one selected and the hit point. The second rule 
is needed only if this hit point happens to be at the bottom of a non-empty column, and 
merging is need: merging between two consecutive columns being possible only when the 
entries of the column on the right are all smaller than the entries of the one on the left. 
Figure 9 illustrates all covers, of a given parking function, for this order. 




For n = 3 and r 
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Figure 9. Covers of 35661100 (for r = 



1). 
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Lemma 3. If ' ip < ip' in r-Tamari order, for r-parking functions ip and (p' , then st(p) = 
st (<//) . The smallest such parking functions is obtained by filling the shape (8n + l n )/Sn^ 
with the permutation a = st(<£>) _1 . 

(r) 

In other words, each of the sets T> a is equipped with an order, which coincides with the 
Tamari order when o is the identity permutation. 

The (/-inversion statistic. Let us recall, from [11], the notion of "dinv" statistic on 
parking functions, which consists in counting their "diagonal inversions". For b > a, 
consider two entries a = r(i,j) and b = r(i',j'), in such a standard tableau r, of shape 
(A + 1")/A, respectively sitting in the cells (i,j) and (i',j')- The pair (6, a) is said to be 
a diagonal inversion if either 

(i) a and b lie on the same diagonal, with a to the right of b. Equivalently, this to 
say that 

i' + f — i + j and i > i' , 

(ii) or, b lies on the "next" diagonal from that of a, with a to the left of b. Equivalently, 
this to say that 

i' + j' = (i + j)+l and i < i . 

In Figure 8, the pairs (8,2), (8,4), (6, 1), (6,3), (7, 1), and (7,3) are diagonal inversions 
of the first type; whereas the pairs (4, 1) and (2, 1) are diagonal inversions of the second 
type. Observe that this notion of diagonal inversion also makes sense for semi-standard 
tableaux of shape (A + l ra )/A. 

We now extend the notion of diagonal-inversion to the context of r-parking functions. 
We first observe that, for a fixed r, two lxl cells (in NxN) intersect several lines of 
slope — 1/r. Indeed, consider two such cells with south-west corner sitting respectively 
in position (i, j) and (i',j')- The number of consecutive diagonals, with equation of the 
form x + r y = k (fc 6 N), that cross the "west wall" of the two cells, is equal to 

max(0, r — \k — k'\), 

where k = i + rj and k' = %' + rj'. In this expression, r is the number of such diagonals 
that cross the west wall of the cell containing 6; whereas (in the case where it makes 
sense) \k — k'\ is the number of such diagonals that only cross the west wall of the cell 
containing b. Hence the difference counts the right number of diagonals, as illustrated 
in Figure 10. 

The r-diagonal inversions occur on such diagonals, as well as on "near" diagonals. Tech- 
nically, consider a r-Dyck path as a partition A contained in the r-staircase Sn ■ For 
a < b, let a = r(i,j) and b = r(i',f) appear as entries in a (semi-) standard tableau r of 
shape (A + l n )/A. The pair (b, a) counts as a r-diagonal inversion, as many times as it 
occurs along a line x + r y = k crossing both west-walls of the cells in which they appear, 
when i > i'; and as many times as a similar property occurs for near diagonals, when 
i < %' . In formula, set k = i + rj and k! = i' + rj', then we have 
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Figure 10. Diagonals of slope — 1/r crossing two cells of (A + l n )/A. 

(i) max(0,r — \k' — k\) r-diagonal inversions, if i > i', and 

(ii) max(0,r — \k' — k — 1|) r-diagonal inversions, if i < i'. 

To extend this to the semi- standard case, we consider that a = b contributes for as many 
r-diagonal inversions, as the minimum of the two values 

max(0, r — \k' — k\) and max(0, r — \k' — k — 1|). 

We denote by dinv r (r) the number of r-diagonal inversions of a semi-standard 
tableau r. 

In the correspondence between standard tableaux and parking function, let us replace 
standard tableaux by semi-standard tableaux. The result is said to be a r-semi-parking 

(r) 

function, and we denote by WA and the set r-semi-parking functions of height n. 
Similarly, we denote by W a the set of r-semi-parking functions that are of shape a, for a 
given r-Dyck path a. Observe that both these sets are infinite, since there is no bound 
on the possible values of the entries. 



Diagonal descents. We now aim at understanding how certain inequalities between en- 
tries of a semi-parking function are sufficient to characterize all positions where diagonal- 
inversions occur in this tableau. This is the reason behind the introduction of the "diago- 
nal descent set" of a parking function. Consider the reverse r-diagonal lexicographic 
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Figure 11. Diagonal reading permutation 82476135, for r = 2. 
order oeNxN defined by 

iff 



<d (i',f) 



% + rj > i' + rj' 



or, 



i + rj = i' + rj', and i < i'. 
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Reading the entries of a semi-parking function ip in increasing <d order, we obtain its 
r-diagonal reading word w r (ip). Informally, we read each diagonal from left to right, 
starting with the one that is furthest from the origin. Clearly, when ip is a parking 
function, this reading word is a permutation of n. For example, 82476135 is the r- 
diagonal reading permutation of the parking function of Figure 11. 

The r-diagonal descent set of a r-parking function (p, denoted by Desc r (<£>), is the set 
of entries a of <p, for which a + 1 appears in a higher position than a in <^ order. In 
other words, a + 1 appears to the right of a in the r-diagonal reading word of ip. For 
example, {1, 3, 5, 6, 7} is the descent set of the 2-parking function of Figure 11. 

Recall that the descent set of a permutation a is 

D(a) := {i |1 < i < n, a { > cr i+1 }. 

The r-diagonal descent set of a semi-parking function ip may be described in terms of 
this classical notion as: 

Desc r (y?) = D(st(w r ((p))~ 1 ). 
This is simply the descent set of the inverse of st(w r (<p)). 

Array encoding. To better manipulate all these notions, it is useful to encode height 
n semi-parking functions as 2 x n arrays, in the following manner. Let us label diagonals 
of the augmented r-staircase shape 5n^ + l n , by their distance from the outermost r- 
diagonal. Thus a cell belongs to the fc th -diagonal if and only if 

r(n — 1 — j) — i = k. 

Such a labeling is illustrated in Figure 12. 
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Figure 12. Diagonal labeling, for n = 5 and r = 3. 

We construct a vector (t>i, t> 2 , . . . , v n ) by listing the entries of a semi-parking function, 
reading down columns starting with the leftmost, and going from left to right. In the 
same order, we record in a word (ui, m 2 , . . . , u n ) the diagonals to which each entry belongs. 
The array encoding of a semi-parking function (p then consists in 

vi v 2 ■■■ v n 

ip = 

\Ui u 2 ■ • ■ u n _ 

Since we can clearly recuperate its original description from this data, we use the same 
notation both for the semi-parking function and its array encoding. Observe that the r- 
semi-parking function condition translates to the following two constraints on this array. 
For 1 < i < n, we have 
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(i) < Ui < Ui-i + r, and 

(ii) if Ui = + r then Vi > Vi-i. 

For example, the encoding that corresponds to the 2-parking function of Figure 11, is 

8 6 7 5 1 2 4 3" 
02132002' 

Two consecutive entries and sit in the same column, exactly when Ui = + r. 
This shows why we must have condition (ii) above. Also, the r-diagonal reading word 
w r ((p) of ip is easily obtained by successive left to right readings of the top row, according 
to increasing values of the bottom one. In our running example, we first read off 8, 2, 4 
that sit above 0's; then 7, sitting above a 1; then 6,3, sitting above the 2's; and finally 
5 which sits above 3. Our previous area statistic is simply the sum of the uf. 

n 

area(y?) = U; L . 
i=i 

5. Polya/Frobenius enumeration 

It is interesting to take into account the action of S n on r-parking function of height n. 
The typical tool for this is the Polya cycle index series, or equivalently the Frobenius 
characteristic. The first notion is usually considered in the context of permutation ac- 
tions, but it can be viewed as a special case of the Frobenius transform of the character 
a linear representation of the symmetric group. To ease future considerations, we will 
adopt this second more general approach. 

Symmetric functions. To simplify the presentation, we assume that w = W\,W2, w^, . . . 
is a denumerable set of variables. The ring A of symmetric polynomials in the variables 
w is graded by degree: 

A = 0A, 

d 

where denotes the degree d homogeneous component. It is easy to see that A^ affords 
as basis the set {m^w)}^ of monomial symmetric polynomials. Recall that these are 
indexed by partitions A = A1A2 • • • A& of d. Then 

™a(w) = w h w i2 ■ ■ ■ 

ill •••)** 

where the sum is over all possible choices of k distinct indices. For example, we have 
m 2 i(x, y, z, . . .) = x 2 y + x 2 z + y 2 z + x y 2 + x z 2 + y z 2 + . . . 



'We will often say functions to underline that we work with "polynomials" in infinitely many variables. 
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We mostly use the notation of [21], so that /ifc(w), efc(w), pfc(w) respectively denote 
the complete homogeneous, elementary, and power sum symmetric polynomials. 
These are characterized by the following generating function identities 

fc>0 i 

^e fc (w)t fc = J (1 + and (27) 



fc>0 



E^ = "EMI-™.*). (28) 

fe>l i 

Recall that, setting 

/i A (w) := /iAx(w) • • -h\ k (w), e A (w) := e Al (w) • • -e Afc (w), and 
p x (w) := PAi(w)---pA fc (w), 

the three sets {h\(w)}\i- d , {e\(w)}x\-d, and {PA( w )}Ahd are bases of A d . The classical 
Cauchy kernel identity (see Macdonald's book [21]) states that 

hi a» 

holds, if and only if {/^(w)}^ and {^(w)}^ form a dual pair of bases for the Hall 
scalar product, that is 



C/a(w),^(w)) 



1 if A = fj,, 
otherwise. 



The two basis {pa( w )}a and {p A (w)/^ A } A constitute such a dual pair, if z\ denotes the 
integer 

z x := l dl d^ 2 d 2 \---n d " dj, 

with di = di(X) equal to the number of parts of size i in A. The monomial and complete 
homogeneous basis are dual to one another. The Schur functions s A (w) are characterized 
by the fact that they form an orthonormal basis 



(s A (w),s M (w)) 

and that they expand as 



1 if A = /i, 
otherwise, 



s A (w) = m A (w) + 22 K \v m A*( w )> (3°) 

with coefficients known as the Kostka numbers. Recall that these are the positive 
integers obtained by counting the number of semi-standard tableaux of shape A, and 
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content fi. This is to say, tableaux that are filled with fii copies of the integer i. In 
particular, the partition /x has to be smaller than A in dominance order, so that 



Indeed, if j > i, the integer i cannot appear in row j; thus all copies of integers smaller 
or equal to i must placed in the first i rows. For example, we have 



S 4 (w) = 7714 (w) + TO 3 l(w) + m 22 (w) + 771211 (w) + 777mi(w) 

S3i(w) = m 3 i (w) + m 22 (w) + 2 m 2 n(w) + 3 m 1 m(w) 
«22(w) = m 22 (w) + m 2 n(w) + 2m im (w) 
S2ii(w) = 777211 (w) + 3 rami (w) 
smi(w) = rami(w) 



Observe that s n (w) = h n (w), and Sn-i( w ) = e n (w). As usual uj stands for the in- 
volution that sends the Schur function s M (w) to the Schur function s M '(w), indexed by 
the conjugate partition. Thus uh n (w) = e„(w). It may also be worth recalling that 
wft(w) = (-l) j -Vj(w). 



Frobenius characteristic. Given a linear action of §„ on a rf-dimensional vector space 6 
V, the associate Frobenius characteristic is defined as 



where M a is the d x d matrix corresponding to the action of a on V, and A(cr) is the 
partition corresponding' to the cycle structure of a. Observe that, when S n acts by 
permutation on a set E (letting V := CE be the corresponding free vector space), then 
the trace of the matrix M a is simply the number of fixed points of a acting on E. Observe 
also that dim(V) appears in V(w) as the coefficient of p™/n\. 

The Frobenius characteristic establishes a correspondence between linear representations 
of S n , and symmetric functions. This correspondence is linear, and sends irreducible 



'In our case, this could be a free vector space, on some set E on which S n is acting by permutation. 
Its parts are the cycle lengths. 



Hi > Ai, 

Mi + \H > Ai + A 2 , 

Mi + M2 + M3 > Ai + A 2 + A 3 

etc. 




(31) 
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representations to Schur functions. More specifically, we have 

(Vi © V 2 )(w) = Vi(w) + V 2 (w), (32) 

(Vi <g> V 2 )(w) = V 2 (w) ■ V 2 (w), (33) 

V(w) = s A (w) iff V is irreducible, (34) 

V(w) = s n (w) iff V is the trivial representation, (35) 

V(w) = sn...x(w) iff V is the sign representation, (36) 

dimV = i<V(w),tf). (37) 

Recall that there are exactly as many irreducible representations of § n , as there are 
conjugacy classes in S n , and that these are naturally in bijections with partitions. Indeed, 
a conjugacy class is entirely characterized by the cycle structure of permutations that lie 
in it. 

Let {V\}\\- n be a complete set of representatives of irreducible representations of S n . The 
fundamental theorem of representation theory, stating that any representation decom- 
poses uniquely into irreductibles, may thus be casted equivalently as 

V = 0a A V A iff V(w) = ^a A s A (w), (38) 

Ahn Ahn 

where the positive integers a A give the multiplicity of the irreducible representation V A in 
V. In other words, the symmetric function V(w) contains all the necessary information 
about the representation V: its dimension, how it decomposes into irreducible, etc; in 
an easy to retrieve fashion. 

Recall that dim(V A ) is the number of standard tableaux of shape A, given by the well 
known hook formula 

A = || j < 39 > 

ll(ij)eA n v 

where we attach to each cell in A, its hook length 

h i:j := Aj+i + A' m — j — 

Thus, when we have (38), 

dim(V) = ^a A / A . 

Ahn 

Exploiting the fact that the trace of conjugate matrices are equal, we may rewrite (31) 
as 

V(w) = J>(aO^, (40) 

/ihn ^ 

where xif 1 ) stands for Trace(M CT ), for any permutation a of cycle type \i\ and where 



z 



:= l dl d l \2 d ^d 2 \---n dn d n \ 1 



if /i has di parts of size i. 
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An example. To illustrate this notion of Frobenius characteristic, let us consider the 
ring TZ n : = C[x] of polynomials in the n variables x = xi,X2, . . . ,x n ; which affords the 
linear basis of monomials 

Y <» ._ ~ai„a-2 . . . ~On 
- — x l x 2 x n i 

indexed by n-vectors a = (ai, 0,2, ... , a n ) in N n . This ring is graded with respect to the 
degree 

deg(x a ) := |a| 

:= a% + a 2 + . . . + a n , 

so that we have the decomposition 

where TZn stands for the homogeneous component of degree d of TZ n - In other words, 
TZ^ is the span of the monomials of degree d. The Hilbert series of TZ n records the 
dimension of each of these component in the form of the power series 

d>0 

It is no hard to show directly that 

K n (q) = (t~~) ' or equivalently dim(^) = (™ + * ~ 1 \ (41) 

However, we will actually prove a more general formula. 

The symmetric group acts on lZ n by permutation of the variables: 

Although lZ n is infinite dimensional, each of the component TZn is finite dimensional. 
We exploit this to define its graded Frobenius characteristic 

^ n (w;g):=^^(g)g d . (42) 

d>0 

This makes sense, since the action of S n respects degree. In fact, this is a permutation 
action on the basis of monomials. From the definition, we may calculate an explicit 
formula for TZ n (w;q) as follows. A straightforward reformulation of (42) gives 



Kn(w;g) = £ E ( 43 ) 

fihn cr M -x a =x a " 



with (Tp any fixed permutation of cycle type fi. We thus have to count monomials, such 
that cr^ • x a = x a , with weight g' a ' . To do this, we select to be 

cr M = (1, . . ., /ii) {tlx + 1, . . ., ill + /i 2 ) • • • (/ii + ... + /ifc-i + 1, . . ., Ill + . . . + /i fc ), 



22 



F. BERGERON 



if fi = //i/i2 • • • l^k- Clearly, er M ■ x a = x a if and only if 

ai = --- = a w (=^i), 
GVi+i = • • • = a w + M2 (= 62), 



o Ail +...+ Atfc _ 1 +i = • • • = fl Ml +...+ w (= &fc)- 
The coefficient of p fJi (w)/z fl in (43) is thus seen to be equal to 

k 



and we get the final formula 

M hn ^ i=l 1 9 



(6i,...,6 fc )eN* 1=1 



fc 



The resulting expressions may be expanded in terms of Schur functions, and we calculate 
that 

Ki(w;q) = 1 si(w), 
(!-?) 

^(w;g) = ri- (s 2 (w) + gs u (w)) , 

(1 - g)(l - g 2 ) 

^s(w;g) = _ _ 1 g 2)( 1 _ g3 ) M w ) + ( g + g 2 ) g 2i(w) + q 3 sm(w)) , 
^(w;g) = (1 - - - ^ _ qA) ( S4 (w) + (g + g 2 + g 3 ) S3 i(w) 
+ (g 2 + g 4 )s 22 (w) + (g 3 + g 4 + g 5 )s 211 (w) + g 6 s llu (w)). 

The nice numerators of these expressions are in fact special instances of the well known 
Hall-Littlewood polynomials. The coefficients of their expansion in the Schur basis 
are known as the g-Kostka polynomials. The fact that these lie in N[q] was explained 
combinatorially by Lascoux and Schtzenberger (see [20]), in terms of the "charge statis- 
tic" on standard tableaux. 

Parking function representation. For a fixed integer partition A, consider the tran- 
sitive permutation action of § n on the set of A-partitions of {1,2, . . . ,n}, i.e.: having 
parts size specified by A. For example, with A = 32, the corresponding set contains the 
10 partitions 

{{1,2,3},{4,5}}, {{1,2,4},{3,5}}, {{1, 2, 5},{3,4}}, {{1, 3, 4},{2, 5}}, 
{{1,3,5}.{2,4}}, {{1,4,5}.{2,3}}, {{2, 3, 4}.{1, 5}}, {{2, 3, 5},{1, 4}}, 
{{2,4,5}.{1,3}}, {{3,4,5}.{1,2}}. 
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Then, one checks that the associated frobenius characteristic is h\(w). This is can be 
done by a straightforward calculation; but one may also first observe that (35) states 
that this holds for n-partitions, and then judiciously apply (33) to get the general case. 
The dimension of this representation is clearly given by the multinomial coefficient 



n \ n\ 



\J Ai!A 2 !---A fc ! 

For a given r-Dyck path a of height n, it easy to see that the the action of S n , on the set 
of parking functions of shape a, is isomporphic to the action of S n on A-partitions. Here 
A is the partition obtained by sorting the parts of the composition 7(a) in decreasing 
order. Modulo some calculations, we have (see [24]). 

(r) 

Proposition 1. The Frobenius characteristic of the space Vn , spanned by r-parking 
functions of height n, my expressed in the following equivalent manners 



VP(w) = *7(«)( w ) ( 45 ) 

1- V(m + 1)^) ^ (46) 



rn + 1 f-^ z x 

Ahn 



^n(" ! ! A, V^(w) (47) 



rn + 1 ' \ n 

Ahn 



^TT 5>(l rn+1 ) *a(w) (48) 



rn + _ 

Ahn 



Ern (rn — 1) ■ ■ ■ (rn — £(X) + 2) , . , , . 



Ahn 



where di(X) is the number of parts of size i in A. 

Recall that s\(l, 1, . . . , 1) may easily be calculated using the formula 

s a( i__i) = n ( 5 °) 

N copies (m)GA V 

where stands for the hook length of the cell in A (here identified with its Ferrers 
diagram). Thus, we may calculate that the respective multiplicities of the trivial and 
the sign representation in Vn \ are given by the Fuss-catalan numbers 

-^Ul"'«)~( ( ' + 1) "), (51) 
rn + 1 rn + 1 \ n J 

and 

e n (l rn+1 ) = 7 ^-r( rn )- (52) 



rn + 1 (r — l)n + 1 \ n 

In the sequel we will see that a "sign-twisted" version of the Frobenius characteristic 
Vn r (w) plays an important role in the study of "diagonal harmonics". Technically, this 
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sign-twisting simply means that we replace ^ 7 ( a )(w) by e 7 ( Q )(w) in (45). In fact, a 
weighted enumeration of semi-standard tableaux of shape (A + l n )/A is of interest in 
that context, and we will be considering the expression 

V n (w; q, t):=J2Yl V inv Mw T , (53) 

XCS n T 

where the second sum runs over semi-standard tableaux r, of shape (A + 1")/A, and w T 
is the tableau monomial 




(ij)6(A+l»)/A 

in the variables w. Observe that |5 n /A|, which is the number of cells in the skew shape 
5 n /X, is simply the area of the Dyck path corresponding to A. 

We tie this with the current section by observing, as discussed in [14], that we have 

J^w T = e 7(A) (w). (54) 

T 

Hence, 

V n {w;q,l) = ^/"/ A le 7(A) (w). (55) 

AC5 n 

Along these lines, let us also mention that Theorem 3.1.3 of [14] states that 

V x n (w;t):= ^^w., (56) 

t:(A + 1"-)/A >N 

t semi-standard 

is symmetric in the variables w, and that it is Schur positive. This is to say that 
V*(w;t) lies in N[*]{s A (w) | A h n}. 



Parking and interval representation. Just as for parking functions, the symmetric 
group acts on pairs in the set Qn (see (22)). We have the following. 

Proposition 2. The Frobenius characteristic of the space Qn^ (spanned by pairs (a,tp), 
with <p a r-parking functions of height n, and a < i((p) in r-Tamari order), may be 
expressed in the following equivalent manners 

EfiSftw) = E E X(a < P) h lW) (w) (57) 

= (^Tiy exp l v E (rn+1) l k )—) (58) 

1 v^tt rn+l f(r + l)(rn + j + 1)\ 

K * Ahn j'GA J 7 
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Proof. Equality between (57) and (58) is shown in [8]. To check that (58) equals (59), 
we start with the generating function identity (established below) 

\fc>i ^ ' J j>i J ' 

where we use the (shifted) Pochhammer symbol notation 

(u) (j) := (u-1) ( u -2)---(u-j + l). 

In light of (26) and (28), one way of parsing Formula (60) is to consider it as stating 
that hj(y) evaluates to ab (ab + (j— i) , when Pk(y) is replaced by a( k )- On the other 
hand, using Cauchy identity (29), we have that 

exp ( J^Pfc(y) ) = hx ^ m ^(w). 

\fc>l / A 

Hence, setting a = rn + 1 and b = r + 1, we get 



exp Z> n + 1) v k — 
\k>i v 7 



Y[(r + 1) (rn + 1) ((r + l)(rn + j + l)) (j) m A (w). 
a jex 

Since it is straightforward that 

rn + 1 fir + 1) (rn + j + 1) 



(r + 1) (rn + 1) ((r + l)(rn + j + 1)) 



U) 



rn + j + 1 \ j 
we have finished our proof. □ 

To establish (60), we proceed as follows. For any positive integer b, consider the series 
(obtained by setting a = 1 in the left-hand side of (60)) 



We first show that Z(i) satisfies the equation 

Z(t) = (l+tZ(t)) b . (61) 
Equivalently, the logarithmic derivative of Z(t) is such that 

Z'(t) bZ(t) 



Z(t) l-(b-l)tZ(t) 
and we want to check that this is equal to 



(62) 



no-Eft)*- 

fe>l v 7 
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Now, it happens that Y(t) satisfies the equation 

Y(t) - (Hl + tY(t))) b m 

Replacing Y(t) by the right-hand side of (62) in this last equality, after calculations, we 
see that Z(t) satisfies (61). Using Lagrange inversion, it is then easy to check that the 
right-hand side of (60) is indeed equal to Z(t) a . 

6. Diagonal harmonics 
In the ring TZ n : = C [X] of polynomials in the three sets of variables 

X := 




"K n is the subspace of polynomial zeros of the polarized power sums differential op- 
erators 



n 



P a (dX) :=^dx a 3 dy b 3 dz c p (64) 

3=1 

for a = (a, b, c) running thought the set of vectors of norm |a| := a + b + c, with 1 < 
|a| < n. Here, dv denotes derivation with respect to v. It is clear from the definition 
that "Kn is closed under derivation. For one set of variables, say x = Xi, . . . , x n , 0~C n is 
entirely described by classical theorems (see [18, section 3.6]) as the span of all partial 
derivatives of the Vandermonde determinant 

A n (x) = det (x?). . 

y ' \ 1 / l<i<n, 0<j<n— 1 ' 

It is of dimension n\, and isomorphic to the regular representation of S„, for the action 
that permutes variables. The description of the case of two sets of variables took close to 
15 years to be finalized (see [16]), and much remains to be understood. It is of dimension 
(n + l) n_1 . It has been observed by Haiman (see [15]) that the dimension of "K n seems 
to be given by the formula 

dim^{ n = 2 n (n + l) n " 2 , (65) 
but, until very recently, almost no one had further studied the trivariate case. We have 
endeavored to do so, considering more general spaces (see definition (80)) for which 
we have experimentally found out that 

dimftW = (r + l) n (rn + l) n " 2 . (66) 

Symmetric group action. We turn both TZ n and "K n into § n -module by considering 
the diagonal action of symmetric group § n on variables. Recall that, for a G §at, the 
polynomial a ■ f(X) is obtained by replacing the variables in f(X) by 

( X a(X) X a (2) ■ ■ ■ 
2/(7(1) y<r(2) ■ ■ ■ 

K z a{ i) ' ' ' 
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We denote by X A the monomial 

X A := {xT---x^).{ V ^...y h -).{zf 
for a matrix of integers 




A = 61 6. • • • 6. 



71 



Writing Aj for the j th -column of A, we may define the degree of a monomial X A to be 
the vector 

71 

deg(X A ) :=5>„ 
j=i 

in N 3 . Clearly the diagonal action preserves degree. The total degree tdeg(X A ) of a 
monomial X A is the sum of the components of deg(X A ). An interesting subspace of 0~C n 
is the space A n of alternating polynomials in "K n . Recall that these are the polynomials 
such that a ■ f(X) = sign(a) f(X). It is also worth recalling that, in the case of two 
sets of variables, this subspace has dimension equal to the Catalan number C n . Another 
observation of Haiman (loc. sit.) is that, in the trivariate case, we seem to have 

dim.A w = -j-— ( 4ri + 1 ), (67) 
n (n + 1) \ n — 1 J 

More generally, our experimental calculations suggest that the alternating component 
An of has dimension 

dim^'= ' r+1 ' f (r + 1) '" + r V (68) 
n(rn + l)\ n-l J v ' 

Degree grading. The space !K n is graduated by degree, and its homogeneous compo- 
nents CK n)( j (with d G N 3 ) are S^-invariant, hence A n is also graded. There corresponds 
direct sum decompositions 

M n = K nA , and A n = A nA . (69) 

A striking fact is that the space "K n is finite dimensional. In fact, the homogeneous 
component Ji n ,d is non vanishing only if 

|d| = d\ + d 2 + d 3 < 

where d = (di,d2,d 3 ). We thus have the Hilbert series (polynomials): 

(;) (3 

^n(q) := dim (Mn,d) q d , and A n (q) := £ dim(A n , d ) q d (70) 

|d|=0 |d|=0 

of "K n and A n , writing q d for q^qpqf* ' ■ The Hilbert series 3~C n (q) is a symmetric poly- 
nomial in qi, q 2 and q 3 . In fact, using representation theory of GL 3 , one can see that it 
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is always Schur positive. This is to say that it expands as a positive integer coefficient 
linear combination of the Schur polynomials s At (q) in the variables q = qi, q 2 , q 3 - 

For example, we easily calculate that CK 2 affords the linear basis 

{l,x 2 -x 1 , y 2 -yi, z 2 -zi}, (71) 

whose subset {x 2 — x%, y 2 — yi, z 2 — z±} is clearly a basis of A 2 . We thus get the Hilbert 
series 

IK 2 (q) = 1 + (qi + q 2 + q 3 ) = 1 + Si(q), and -A 2 (q) = qi + <? 2 + q 3 = Si(q). 

To calculate ^K„(q) for larger n, we exploit the fact CK n is closed under taking derivatives, 
as well as applying the operators 



i=l 



Here u, v stand for any two of the three sets of n variables: 

x = xi, . . . ,x n , y = 2/1, ...,y n , and z = z lt . . . ,z n . 

It easy to verify that the Vandermonde determinant A n (u) belongs to !K n , whether u be 
x, y of z. Using all this, we can calculate that 

JC 3 (q) = l + 2{ qi + q 2 + q 3 ) 

+2 (ql + ql + ql) + 3 (q x q 2 + q x q 3 + q 2 q 3 ) 

+qf + <?2 + qI + + Q1Q3 + q\q\ + qiql + q\q 3 + q*q\ + q^qs 
by checking directly that we have the following respective bases Ed for each 'K n ^: 
£300 = {A 3 (x)}, S 200 = {^ 1 A 3 (x),^ 2 A 3 (x)}, S 100 = {^A 3 (x) A^AaCx)}, 
03 O3O = {A 3 (y)}, 23 020 = {d yi A 3 (y), ^ 2 A 3 (y)}, a3 010 = {d 2 yi A 3 (y), d yi d m A 3 (y)}, 
S 003 = {A 3 (z)}, S 002 = A 3 (z), ^As(z)}, S001 = A 3 (z), Xl & 2 A 3 (z)}, 
£210 = {£ yx A 3 (x)}, S 120 = {£ yx £ yx A 3 (x)}, S 201 = {£ zx A 3 (x)}, 
£102 = {£ zx £ zx A 3 (x)}, S 021 = {£ zy A 3 (y)}, S i2 = {£ zy £ zy A 3 (y)}, 
Sue = {E yx d Xl A 3 (x),E yx d X2 A 3 (x.),E$A 3 (x)}, 
S101 = { J E zx ^ 1 A 3 (x),E zx 9 :E2 A 3 (x),4 x ) A 3 (x)}, 
Son = { J E zy a 2/1 A 3 (y), J E zy ^ 2 A 3 (y),Eg ) A 3 (y)}, 
S 000 = {1}, Bin = {£ Z x£ yx A 3 (x)} 

Collecting previous observations, and doing some further explicit calculations, we get 
that 

fti(q) = 1, 

ft 2 (q) = l + si(q) 

J£ 3 (q) = l + 2 Sl (q) + 2 S2 (q) + Sll (q) + 3 3 (q), (72) 
J{ 4 (q) = 1 + 3si(q) + 5s 2 (q) + 3s n (q) + 6s 3 (q) + 5s 21 (q) + s m (q) 
+5 s 4 (q) + 4 s 3 i(q) + 3 s 5 (q) + s 4 i(q) + s 6 (q). 
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If we specialize one of the parameters (say g 3 ) to 0, we get back the graded Hilbert series 
of bivariate diagonal harmonics (of overall dimension (n + l)™" 1 ) which has received a lot 
of attention in recent years. In other words, the z-free (or x-free, or y-free) component 
of "Kn coincides with the "usual" space of diagonal harmonics (in the bivariate case, 
as previously considered in the literature). Hence, all of our formulas involving the 
parameters q specialize to known formulas by the simple device of setting one of the 
three parameters in q equal to 0. Evaluating CK n (q) at ^ equal to 1, we clearly get the 
overall dimension of "K n . In other words, the resulting evaluation of (72) agrees with 
Formula (65). 



Graded character. We refine the dimension analysis by taking into account the de- 
composition into irreducibles of the homogeneous components of CK n , see (69). This is 
all encompassed into the graded Frobenius characteristic of !K: 



^)(q;w) :=^q d <i 

deN 3 



n;dl W J 



For reasons already alluded to, the expansion of CK n (q; w) in terms of the S\(w) affords 
Schur positive coefficients in the s At (q). This is to say that 

*4 r) (q; w) = (E<!^(q) ) &(w), n x ,„ e n. (73) 

Ahn \ u / 

Hence there are two kinds of Schur function playing a role here. To emphasize this, we 
denote those in the q- variables by a lower case "s" . For example, we have 

9ti(w;q)=Si(w), 

:K 2 (w; q) = ^(w) + si(q) 5"ii(w), 

^3(w; q) = S 3 (w) + (si(q) + s 2 (q)) S 2 i{w) + (s n (q) + s 3 (q)) S ni (w), 
M 4 (w; q) = S 4 (w) + (si(q) + s 2 (q) + s 3 (q)) S 3 i(w) 
+ (s 2 (q) + «2i(q) + S4(q)) 5"22(w) 

+(*n(q) + s s(q) + s 2 i(q) + s 4 (q) + s 31 (q) + s 5 (q)) S 2 n(w) 
+ (sm(q) + s 31 (q) + s 41 (q) + s 6 (q)) S'm^w). 

The Schur functions S\(w) act as placeholders for irreducible representations of S„. Thus 
the multiplicities of the irreducible representations of § n in Jt n appear as the coefficients 
of the ^(w)^. Although we think that these coefficients should evaluate in the three 
parameters q = (71,(72,93, they are in fact "universal" expressions (see [1]). This is to 
say that, in any k sets of n variables, k = 1, 2, 3, . . ., the above formulas evaluate to the 
"correct" value when we set q = qi, g 2 , . . . , q^. 



(74) 
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Using Formula (50) we can exploit this, for k generic, to get the following polynomial 
expressions in k, when we set all qi — 1. 

M 1 (w;l fc ) = S 1 (w), 

K 2 (w;l k ) = S 2 (w) + kS ll (w), 

:K 3 (w; i*) = 5 3 (w) + (k + ft 1 )) 5 21 (w) + (© + ( fc + 2 )) 5 m (w) 
:K 4 (w; l fc ) = 5 4 (w) + (k + ft 1 ) + ft 2 )) 5 3 i(w), (75) 
+(R 1 )+2(?) + f| 3 )))&(w) 

+(© + m)+2m + m +3m + 
+(©+3( fc : 2 )+4m + ( fc r))^im(w). 

Thus, if we set k = 1, we get the Frobenius characteristic of the space of classical 
harmonics; at k — 2, we get that of the space of bivariate diagonal harmonics; at k — 3, 
trivariate diagonal harmonics; etc. 

Considering only the coefficients of Sii...i(w), we obtain the formulas 

A 1 (l k ) = 1, 
A 2 {l k ) = k, 

Mi k ) = © + m, 
Mi k ) = ©+3( fc : 2 )+4( fc : 3 ) + (T)' 

which respectively give the appropriate values 1, C n , and (67); at k 
k = 3. 



(76) 



1, = 2 and 



Bivariate diagonal harmonics. Before going on with our presentation, let us recall 
some classical facts about the spaces of harmonic polynomials, and recent results about 
the space of bivariate diagonal harmonic polynomials. 

The classical space of harmonic polynomials is n\ dimensional, with graded Frobenius 
given by the formula 8 

n 

X n (w;q) = J(l-q k ) ^ s A (l, g, g 2 , . . .) 5 A (w), (77) 



k=l Ahn 
n 



= na-«*)Enf^- < 78 > 

k= - 

It follows that its Hilbert series is 



Q 

k=l Ahn i 



k=i - - 9 

Observe also that when q goes to 1, we get 

IK n (w; 1) = h\. 

8 This is the numerator of lZ n (yv; q) in Formula (44)). 
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As already mentioned, the space of bivariate diagonal harmonic polynomials has dimen- 
sion equal to (n + l) n_1 . Its (bi-)graded Frobenius characteristic 3i n (w;q, t) affords a 
description in terms of an operator which affords "Macdonald Polynomials" as eigenfunc- 
tions. Rather than recall the necessary notions, so that we could present this description, 
we recall from [14] a proposition giving a combinatorial description for < K n (w; q, 1). Ob- 
serve that we are here setting t = 1. A combinatorial description of "K n (w;q, t) in full 
generality is still but a conjecture. We will come back to this in later sections. The 
notion of "higher" diagonal harmonics is explained in the next subsection. 

Proposition 3 (See [14], Prop. 4.1.1, page 10). The q-graded Frobenius characteristic 
of the space of bivariate (higher) diagonal harmonics is given by the formula 

3CST>(w;g,l)= 9 arca(a) e 7(Q) (w). (79) 

This is the "sign twisted" version of the graded Frobenius for parking function represen- 
tation, expressed in Formula (55). We will see that a similar statement seems to hold in 
the trivariate case. 



Higher harmonics. All of these considerations generalize to the spaces of higher di- 
agonal harmonics 

JK< r > := (-4IT 1 ) n (Al-'Xn) 1 , (80) 

where A n (resp. X n ) stands for the ideal generated by alternating polynomials (resp. 
invariant) polynomials in TZ, and the orthogonal complement is with respect to the 
§ n -invariant scalar product defined by 

(X\X°) ' iA = B> (81) 

1 otherwise. 

Here A\ stands for the product of factorials: 

A\:=(a 1 \---a n V)-(b 1 \---bJ)-(c 1 \---c v }). 

On !Kn the diagonal action is twisted by e r ~ 1 , where e stands for the sign representa- 
tion. For n — 2, the space Ji 2 decomposes as a direct sum of its isotypic invariant 

(r) (r) 

component, J 2 , and alternant component, A 2 ■ 

jfW =J« ©4 r >, (82) 

which afford the respective bases 

S^ r) = {{x 2 - x 1 ) a {y 2 - yi )\z 2 - z x ) c \ a + b + c = r-l}, and 
■BA 2 r) = {{x 2 -x l ) a {y 2 -y l )\z 2 -z l y\a + b + c = r}. 
In fact, we have a universal formula, with both k (where q = qi, q 2 , . . . , q^) and r generic: 

^4 r) ( w ; q) = Sr-i(q) ^(w) + s r (q) S'ii(w), 
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From this we get the universal Hilbert series 

^2 r) (q) = s r-i(q) + s r (q), and A$'(q) = s r (q). 
In particular, setting all g$ — 1, we have the dimension formulas 

^d'H^^f^- 2 ), and ^(iW*^" 1 



as polynomials in the parameters k and r. For the trivariate case, this gives 

dim^ r) = (r + l) 2 , and dim.A 2 r) = ( r + 1 K r + 2 ) , 

which agree with Formulas (66) and (68). The above formulas are both encompassed in 
the following conjecture. 

Conjecture 1. The graded Frobenius characteristic ofKn^ is given by the formula 

IKM(w;g,l,l)= x(a<P)q di ^ ] e im (w). (83) 

Observe that, at q = 1, this is a sign-twisted version of (57). For example, at r = 1, 
calculating as in (75) we obtain the following formulas that include as special cases 
both (79) and (83). 

Hi(w; g ,l fc ) = ei(w), 

JC 2 (w;g,l fc ) = e ?(w) + (g+(A;-l))e 2 (w), 

H 3 (w;q,l k ) = e?(w) + (g 2 + (A; + l)g + (A; + 4)(A ; -l)/2)e 2 (w)e 1 (w) 

+ (q 3 + (k-l)q 2 + (k + 2)(k - 1) q/2 + (k — l)(k 2 + Ak - 6)/6) e 3 (w). 

7. Combinatorial conjectures 

The climax of our story concerns explicit (but still conjectural) combinatorial formulas 
for the graded Frobenius characteristic of the spaces of bivariate and trivariate higher 
harmonic. We first recall the following (which involves a r- refined version of (53)). 

Conjecture 2 (Shuffle Conjecture, see [14]). The graded Frobenius characteristic of the 
space of bivariate diagonal higher harmonics is given by the formula 

K^(w,q,t)= Yl ,frca(A) t dinv r (T) w ^ ^ 



Ac5 (r) r:(A+l")/A— >N 

— n t semi-standard 



where area(A) stands for the number of cells of the skew shape 5n^ /A. 

Although we know that the left-hand side is symmetric in q and t; at the moment of 
this writing, there is no direct explanation of this fact for the right-hand side from the 
combinatorial point of view. On the other hand, by definition we know that "Kn \w, q, t) 
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expands positively (with coefficients in N[q, t})) in terms of Schur functions of the vari- 
ables w. The fact that this is also the case for the right-hand side of (84) follows from 
(56). 

From our previous discussion, we may recast the above formula in terms of r-semi-parking 
functions. In this terminology, the conjectured identity (84) becomes 

J£< r >(w; g ,f)= £ iT^^S. (85) 

where all notions are directly inherited through the translation from semi-standard path 
terminology. Our reason for all this is the ease the rest of our presentation. 

Trivariate extension. The most striking aspect of the new developments presented 
here, is that a similar conjecture seems to hold for the trivariate case. The technical 
statement is as follows. 

Conjecture 3 (see [5]). The graded Frobenius characteristic of the space of trivariate 
diagonal higher harmonics is given by the formula 

WfJ(w; q,t, 1) = Y, E *( a ^ ? d(a '^ dmVrM w,. (86) 

Observe that this conjecture is lacking a description of the third statistic that should 
appear here, rather than be specialized to 1. 

Quasi-symmetric functions. Both of the above conjectures involve infinite sums. 
These may be turned into finite expressions using the theory of "quasi-symmetric func- 
tions" . We start with the observation that, for a fixed shape (A + l ra ) / A, all semi-standard 
tableaux of some specific families have the same diagonal-inversion statistic. More pre- 
cisely, the value of dinv r (r) is entirely characterized once some inequalities are known to 
hold between pairs of entries of r, irrespective of the actual values of these entries. We 
will group together all such terms to obtain quasi-symmetric functions. Both conjectures 
then take the form of finite sums of expressions involving quasi-symmetric functions, in- 
dexed by compositions associated to diagonal descent sets. 

Recall that degree d homogeneous component, of the ring Q of quasi-symmetric function, 
is spanned by the monomial quasi-symmetric functions M 7 (w), which are indexed 
by compositions 7 = (ci, . . . , C&) of d: 

M 7 (w) := < 

where we write u Ct w to say that u consists in a choice of k of the variables in w, these 
variables being ordered as they appear in w. In other words, 

u 7 = w^w% ■ ■ ■ w%*, with z'x < z 2 < . . . < 4. 
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For example, 

M 2 i(x, y, z) = x 2 y + x 2 z + y 2 z, and M 12 (x, y, z) = xy 2 + xz 2 + yz 2 . 

It easy to see that the monomial symmetric function m\ is equal to the sum of all M 7 , 
with 7 running through the set of all compositions that give A when their parts are 
reordered in decreasing order. 

The so-called fundamental basis (or F-basis) of Q is defined as the set of quasi- 
symmetric functions 

F 7 (w) :=^M a (w), 

where the sum is over all compositions a that are finer than 7. Recall that this is to say 
that 7 may be obtained from a by adding together consecutive parts. Thus, we have 

#5121 (w) = M 312 i(w) + M 3im (w) + M 2 ii2i(w) + M21HH 

+Mi212l(w) + Mi21Hl(w) + M U 112l(w) + Miiinii(w). 

An alternate description exploits the classical bijection between partitions of d and sub- 
sets of {1,2, ...,d— 1}. Recall that this is the bijection that sends the composition 
7 = (ci, . . . , Cfc) to the set of partial sums 

5(7) := ji> 1 i<j 

Thus we have a ^ 7 in refinement order if and only if S^) C S(a). 

Finite form of conjectures. We can now collect terms in (85) and (86) so that they 
be expressed in terms of the F-basis. The resulting formulas take the form 

X$[\w;q,t) = E ^^(w), (87) 

Mi r) (w;g,U) = Yl J2x(a<P)q d{ ^h d ™^) Fe{ ^ w ). (88) 

with summation now being over the (finite) set of parking functions, and where 0(<p) 
stands for the composition of n encoding the diagonal descent set Desc r (<£>), of the parking 
function ip. 

It is shown in [14] that, for any r-Dyck path (3, the sum 

D fi (w;t) := t d ™ r(v) Fe( v) W 

expands, with coefficients in N[t], in terms of Schur functions in the variables w. More- 
over, at t — 1, we have 

Dp(w] 1) = e 7(/3) (w). 
Thus, Conjecture 3 implies Conjecture 1. 
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From the algebraic point of view, it is clear that we must have the identity 

ftW(w;g,l) = ftM(w; g ,0,l). (89) 

It would thus be very interesting to find a direct combinatorial explanation (independent 
of the conjectures) of the equality that would have to hold between the corresponding 
right-hand sides of (87) and (88). 



8. Other research directions 

Beside those explicitly considered in the text above, several natural questions arise. Chief 
among these is to develop a similar combinatorial setup for other finite reflection groups, 
for which much of the representation theory has been clarified in recent years (see [4]). 

There has also been recently a lot of work on geometric (or polytopal) realizations of 
the associahedron (the Tamari poset for r = 1). This leads to the natural question of 
describing similar constructions for all r-Tamari posets. Figure 6 suggest a tantalizing 
outlook on this. 
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